Thermoelectric spin transfer in textured magnets 
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We study charge and energy transport in a quasi-lD magnetic wire in the presence of magnetic 
textures. The energy flows can be expressed in a fashion similar to charge currents, leading to new 
energy-current induced spin torques. In analogy to charge currents, we can identify two reciprocal 
effects: spin-transfer torque on the magnetic order parameter induced by energy current and the 
Berry-phase gauge field induced energy fiow. In addition, we phenomenologically introduce new 
/3— like viscous coupling between magnetic dynamics and energy current into the LLG equation, 
which originates from spin mistracking of the magnetic order. We conclude that the new viscous 
term should be important for the thermally induced domain wall motion. We study the interplay 
between charge and energy currents and find that many of the effects of texture motion on the 
charge currents can be replicated with respect to energy currents. For example, the moving texture 
can lead to energy flows which is an analogue of the electromotive force in case of charge currents. 
We suggest a realization of cooling effect by magnetic texture dynamics. 

PACS numbers: 72.15.Jf, 75.30. Sg, 72.15.Gd 
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The notion of the Berry phased naturally appears in 
the description of magnetic texture dynamics in the limit 
of strong exchange field The spin up and down with 
respect to local magnetization electrons experience "ficti- 
tious" electromagnetic fields.— These fields have opposite 
signs for spin up and down electrons and result in the 
Lorentz forcej^ It has been realized that the spin-transfer 
torque (STT) is a reciprocal effect to the electromotive 
force (EMF) associated with this Lorentz force.— In real 
systems, the exchange field is finite leading to spin mis- 
alignments with the texture, and more realistic descrip- 
tion should take such effects into account via (3 terms in 
the Landau-Lifshitz-Gilbert (LLG) equation.i^ 

Recently, interest in thermoelectric effects has con- 
siderably increased as new experimental data has been 
available.^ The Peltier effect describes heat transfer ac- 
companying the current flow. The opposite is the See- 
beck effect that describes the thermo-EMF induced by 
temperature gradients. The Peltier and Seebeck thermo- 
electric effects as well as the thermoelectric STTs have 
been studied in multilayered nanostructuresi^ Hatami et 
al. proposed a thermoelectric STT as mechanism for do- 
main wall motion.® Berger and Jen and Berger observed 
and discussed domain wall (DW) motion induced by heat 
currentsji^ Thermal STTs may soon be employed in the 
next generation of nonvolatile data elements for reversal 
of magnetization. Thermoelectric nano-coolers can find 
applications in the nanoelectronic circuits and devices. ii 

In this Rapid Communication, we study continuous 
magnetic systems which can be relevant to DW motion^ 
and spin-textured magnetsj^ We phenomenologically de- 
scribe thermal STTs in a quasi-lD magnetic wire with 
magnetic texture. The Berry-phase gauge field induced 
energy fiow turns out to be reciprocal effect to the ther- 
mal STT and both effects can be formally eliminated 
from the equations of motion by properly redefining the 
thermodynamic variables which is reminiscent of the non- 
dissipative STTsiS- We further generalize our description 




Figure 1: (color online). In quasi-lD magnetic wire, charge 
current density j is induced by potential gradients dx^J,, tem- 
perature gradients dxT and EMF dt^ produced by the Berry 
phase which is acquired by the electron spin following the 
instantaneous magnetic proflle. Coupled viscous processes 
arise once we relax the projection approximation. The mag- 
netic texture ra.{x,t) responds to the effective field H(a;,t). 



by including viscous Hke terms corresponding to spin 
misalignments. These viscous effects turn out to be im- 
portant for the thermally induced DW motion and can 
lead to such effects as cooling by magnetic texture dy- 
namics. We also find that the Peltier and Seebeck effects 
can be modified and tuned by the magnetic texture dy- 
namics. 

Consider a thin quasi-lD magnetic wire with the mag- 
netization m{x,t) (along the spin density) in the presence 
of chemical potential ^{x, t) and temperature T{x, t) gra- 
dients. We would like to construct a phenomenological 
description of our system based on thermodynamic vari- 
ables introduced above and their conjugate forces. The 
ferromagnetic wire is supposed to be thermally isolated, 
and after being perturbed by nonequilibrium chemical 
potential, magnetization and temperature gradients, the 
wire evolves back toward equilibrium according to the 
equations of motion, producing entropy. Note that we al- 
low this equilibrium state to be topologically nontrivial. 
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e.g. a magnetic DW or vortex. The state of partial equi- 
librium can be described by thermodynamic variables Xi 
and their conjugates (generalized forces) Xi — dS/dxi 
with the entropy and its time derivative being: 



i.k—l i—1 



XiXi 



We initially consider the entropy E:{p, pu as a func- 
tion of the density of electron charge p, the density of 
energy pu and the magnetization direction. The mag- 
nitude of the magnetization is not treated as a dynamic 
variable, assuming sufficiently fast spin-flip relaxation. 
The conservation laws of energy and charge provide lin- 
ear relations: 



-dxi , Pu = -dxju , 



(1) 



where we introduced the charge current j and energy cur- 
rent ju ■ For conserved quantities, it is more convenient 
to work with fluxes j and ju instead of densities p and 
Pu which leads to equivalent description due to the linear 
relations in Eqs. HI). 

We will not be concerned with the general expression 
for the entropy but rather concentrate on identifying the 
thermodynamic variables and their conjugates by calcu- 
lating the time derivative of the entropy. Suppose that 
we flx the texture, then the rate of the entropy change 
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(2) 

where we introduced the modified energy current jq = 
ju — fJ'.j that describes the energy flow without the en- 
ergy corresponding to the chemical potential p (from now 
on only the energy current jq is considered). The chem- 
ical potential is defined as a conjugate of the density of 
charge. The introduction of jq is necessary to avoid the 
unphysical gauge dependence of the energy current and 
the associated kinetic coefficients on the potential offset 
for the whole system. We are now ready to write the rate 
of the entropy change for the general case of dynamic spin 
texture: 



§ = (h dx dr 



jq - (f dx ^j - (f dx^- dtm , (3) 
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where in Eq. ([2]) we integrated the term involving jq by 
parts and the conjugate/force corresponding to the mag- 
netization is defined as — 9mS|Q,g = H/T with Q{x) and 
q{x) being the overall charge and energy that passed the 
cross section at point x which corresponds to integrating 
j and jq in time, respectively. As it can be seen from Eq. 
([3]), our other conjugates are —dqS\m,Q = —dx (1/T) and 
— 9Q§|m,(2 = dxp/T. In general, H is not the usual "effec- 
tive field" corresponding to the variation of the Landau 
free-energy functional F[m, p, T] and only when dxT = 
and dxP = the "effective fields" coincide. Let us ini- 
tially assume that even in an out-of-equihbrium situa- 
tion, when dxT ^ and dxP 7^ 0, H depends only on 



the instantaneous texture m(a;). In general, however, we 
may expand H phenomenologically in terms of small dxT 
and dxp. 

In our phenomenological theory, the time derivatives 
of thermodynamic variables are related to the thermo- 
dynamic conjugates via the kinetic coefficients. In order 
to identify the kinetic coefficients, we assume that the 
currents j and jq are determined by the chemical poten- 
tial and temperature gradients as well as the magnetic 
wire dynamics, which exerts fictitious Berry phase gauge 
fieldsfi' on the charge transport along the wire. We then 
have for the charge/energy current gradient expansion: 
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iT^^ + p (^m X dxva + iidxva.^ ■ dtui , 

(4) 



CT^^+p' (mxdxm + p'dxm]-dtm, 



(5) 

where we assume that the coefficients g, ^ and C can in 
general also depend on temperature and texture, for the 
latter, to the leading order, as g = go+r]g{dxm)'^ , etc.. In 
Eqs. ^ and we expand only up to the linear order 
in the nonequilibrium quantities dxP, dxT and dtm and 
to the second order in (3a;m; however, the latter terms are 
expected to be small in practice and only are necessary 
for establishing the positive-definiteness of the response 
matrix. The spin-rotational symmetry of the magnetic 
texture and the inversion symmetry of the wire are also 
assumed to avoid additional and often complicated terms 
in our expressions. Relating dtvci to the generalized force 
—H/T, within the LLG^ phenomenology, we derive the 
modified LLG equation consistent with Eqs. ^ and |[5]), 
with the guidance of the Onsager reciprocity principle: 



s(l + amx)9tm + m x H = —p 



dxva 



/3(m X dx-m) 
■dxT/T, 

(6) 

where we introduced the spin density s so that 
sm = M/7, with M being the magnetization den- 
sity and 7 the gyromagnetic ratio (7 < for elec- 
trons). Equation ^ can be expressed in terms of 
the charge/energy flows by inverting the linear relation 
{j.h} = {9.1,1 h}{dxP.dxT/T}: 



dxp. = -gj + £,jq + p (m X dxvci + pdx-m.) ■ dtra , (7) 



dxT/T = - Oq + (m X dxin + /3'a,m) ■ dtui , (8) 



s(l -I- amx)dtm + m x H = p [c^a^m -I- (3{m x dxm)] j 
+p' [dxm + l3'{mx dxm)] jq , 

(9) 
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where the new coefficients g, ^, C, p, p', (3 and /?' can 
be expressed via the ones marked by tilde, and in Eq. 
([9]), we disregarded the terms of the order ^ {dxrxvYdtm 
contributing to the Gilbert damping. Terms of similar 
order can also appear in case of incompressible charge 
flow and lead to the non-local Gilbert damping.— The 
kinetic coefficients contain information about the con- 
ductivity, a = g, the thermal conductivity, k = l/(Cr), 
and the conventional Seebeck and Peltier coefficients can 
be found from Eqs. (0]) and ^ by assuming j = for 
the former, S = —^/(gT), and by assuming dxT = 
for the latter, 11 = ^/^ — — ^/.g, which also implies that 
g = 1/a + S'^T/k. Equation ^ differs from an ordinary 
LLG equation^ by the extra spin torque terms that ap- 
pear in the presence of the energy flow jq. These new 
torques are similar to the nondissipative and dissipative 
current induced spin torques,— as the former can be re- 
lated to electron spins following the magnetic texture 
and the latter - to electron spins mistracking the texture. 
The phenomenological parameter p (or p = p/ao — p'TIq) 
can be approximated as p/<7o = ph/2e in the strong 
exchange Hmit^ and corresponds to the electron spin- 
charge conversion factor h/2e multiplied by the polar- 
ization p = (ctq — ctq )/(To ) fo = o-Q + ""o ^^"^ ^ minus 
the charge of particles, e.g. for electrons e is positive. 

Similarly, we can consider Eq. ^ under conditions 
of vanishing charge currents and fixed texture, and find 
the spin current resulting from the temperature gradi- 
ents: 2eSsdxT/{l/al + 1/ct^) where Ss = [Sl - S^)/e is 
the spin Seebeck coefficient in the absence of magnetic 
texture. By involving the electron spin-charge conver- 
sion factor again, we can approximate the second spin 
torque parameter p' in Eq. ([9]) in the strong exchange 
limit arriving at 

where we introduced the spin polarization of the Seebeck 
coefficient ps = [S] - S^)/ [S] + S^) = eS,/{2So). When 
the thermal conductivity is mostly due to electron mo- 
tion, we can simplify Eq. lfTO|) further with the help of the 
Wiedemann- Franz Law according to which kq/cto — LT 
where L = 7r^fc^/(3e^) is the Lorenz number. Effects 
such as spin drag^^ can also influence the estimate in Eq. 
(fTO|) . Note that the result in Eq. ifTOl) can also be ob- 
tained from Eq. ([8]) by considering the texture-dynamics 
induced EMF which can lead to the energy currents in 
the absence of charge currents. 

One can redeflne the magnetization (which in turn 
leads to changes in generalized forces dx^/T and 
dxT/T^) in Eq. Q so that the nondissipative parts of 
torques are absorbed into these new deflnitions: 
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where with this choice of thermodynamic variables, Eqs. 
([7]), (HI and ^ will only have spin torque terms propor- 
tional to (3{P'). 

From Eq. ([H, we can write the rate of the entropy 
production: 



dx 
Y 



'2j3pidxm ■ dtm 



'2'Uiq + as{dtmf 
2P'p'jqdxm ■ dtm] 



(11) 



where g ^ go + r]g{dxmy, ^ = + Viidxm)^ and C = 
Co + Vcidxm)'^. Notice that the STTs in Eq. © induced 
by the charge/energy currents can be separated into the 
nondissipative and dissipative parts (/3 terms) based on 
Eq. ifTTj). This separation is, nevertheless, formal as 
in reaHstic metaUic systems the torques will always be 
accompanied by the dissipation due to the flnite thermal 
conductivity k and conductivity a. The dissipation in Eq. 
(fTTI) is guaranteed to be positive-definite if the following 
formal inequalities hold: 



? < y/gC, Vg > 



PpP'p' 
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(12) 

where the first inequality can be rewritten equivalently 
as g > S'^T/k, and should always hold. Other inequal- 
ities are somewhat formal since their proof impHes that 
our theory can describe sufficiently sharp and fast tex- 
ture dynamics for dominating dissipation as opposed to 
the first inequality for proof of which a mere static tex- 
ture assumption is sufficient. Nevertheless, Eqs. (fT2l) can 
serve for estimates of the spin-texture resistivities (?7^ and 
rjg) and the spin-texture Seebeck effect (77^) due to spin 
dephasing. The condition on the spin-texture resistivity^ 
Vi/a- > {PpY / {cisaD also follows from Eqs. (fT2|) . 

Let us now discuss thermal effects that can arise from 
the presence of magnetic texture. In case of a static spin 
texture and stationary charge density, dxj — 0, we can 
write the modification to the Thomson effect by calculat- 
ing the rate of heat generationi^ Q = —dxju from Eqs. 
dH) and jHl): 



= Kd^T + {dTn){dxTf + + T{dTS)jdxT 
+V^dx[{dxmf]dxT + Trjsdx[{dxmf]j , 



(13) 



where k = Ko + rj^ldxm)"^ and S = So + ilsidxm)"^ . Even 
though Eqs. (0]) and ^ are the first order gradient ex- 
pansions, they are sufficient for obtaining Eq. lfT3|) since 
consideration of the second order expansions would only 
lead to even higher order terms in Eq. (fT3|) . The Thom- 
son effect can be modified by the presence of texture and 
some analogue of local cooling may be possible even with- 
out temperature gradients. However, the magnitude of 
the coefficients ?7k and t]s is not accessible at the moment 
and should be extracted from the microscopic calcula- 
tions. 
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Another thermal effect we would like to discuss is re- 
lated to heat flows induced by magnetization dynamics. 
As can be seen from Eq. such heat flows can ap- 

pear even in the absence of temperature gradients and 
charge current flows. When the magnetic texture fol- 
lows a periodic motion, the energy flows should result 
in effective cooling or heating of some regions by specifl- 
cally engineering the magnetic state of the wire and ap- 
plied rf magnetic fields. The spin spring magnets can 
be of relevance.— Alternatively, the texture in our wire 
{i.e. spiral), can result from the Dzyaloshinskii-Moriya 
(DM) interaction^^ relevant for such materials as MnSi, 
(Fe,Co)Si or FeGe. In this case, the end of the spiral can 
be exchange coupled to a homogeneous magnetization of 
a magnetic film subject to rf magnetic field which should 
result in rotation of the magnetization in the film and 
spiral. 

To simulate the spiral rotation and obtain the prelim- 
inary estimates of the effect, we consider the current in- 
duced spiral motion that in turn leads to the energy fiows 
due to the magnetic texture dynamics. The static tex- 
ture corrections due to t]— type terms will be ignored. 
We consider a ferromagnetic wire with the DM interac- 
tion in the absence of the temperature gradients. The 
"effective field" can be found from the Free energyji^ 



J 



(Vm)^ + Fm • (V X m) 



(14) 



where J is the exchange coupling constant and F is the 
strength of the DM interaction. The "effective field" H = 
dmF can be used in the LLG Eq. ([9]). The ground state 
of the Free energy in Eq. lfT4|) is a spiral state m(r) = 
ni cosk-r-|-n2 sink - r where the wavevector k = nsF/ J, 
and rii form the right handed orthonormal vector sets. 
We assume that the wavevector is along the wire, e.g. 
due to anisotropies. As can be found from the LLG Eq. 
([6]), for the case of vanishing temperature gradients, the 
spiral starts to move along the wire in the presence of 
currents and the solution can be described as^^ 

m(r, t) = 171x^+171^ (y cos [k{x — vt)] + z sin [k{x — vt)]) , 

where the x axis points along the wire axi s, rrix = 
{jh/e){l3/a - 1)/(2F - Jk), m± = y/l-ml and v = 
pj{f3/a)sh/{2e). The wavenumber k = 2tt/X and nix 
have been calculated numerically in Ref. |13 and in the 
presence of currents A increases and rrix acquires some 
finite value; however, for an estimate corresponding to 
moderate currents the values given by the static spiral 
ko = F/J should suffice. The maximum current that 
the spiral can sustain without breaking into the chaotic 
motion is jmax ^ 2Te/h:^- 

Using Eq. we are now ready to find the energy 
flow accompanying the current flow as the spiral moves 
with the speed v in the absence of the temperature gra- 
dients: jq = Iloj — m\k'^vp' (3' /(, K, O.SIIoj, where only 
/?'— term contributes to the energy flow in Eq. 
For our estimate, we take parameters corresponding to 



a MnSiii^ the lattice constant a = 0.5 nm, the magne- 
tization density M = A = 20 nm, a = 0.01, 
/3 = /?' = 0.03, cr = 5 X 10^ ^"^'^^^^ = 0.02 eV 
and p — ps — 0.8. By increasing /3(/3') and a and di- 
minishing A and a, the energy flow can be made larger. 
We conclude then that the current-induced magnetic tex- 
ture dynamics can lead to additional energy flows that in 
some cases can be comparable to the energy flows due 
to the Peltier effect. The renormalization of the Peltier 
coefficient should also apply to the Seebeck coefficient 
due to the Onsager reciprocity principle dictating that 
S — n/T. Assuming the absence of the temperature 
gradient, from Eq. (0]), we also flnd correction to the 
conductivity caused by the EMF due to the spiral mo- 
tion: —dxfJ- = j/(To + ni\k'^vpl3h/{2e) wO.Sj/cro. 

Finally, we calculate the speed of the spiral motion in- 
duced by temperature gradients in the absence of charge 
currents. In full analogy to the spiral motion induced by 
the charge currents and using Eq. we flnd the spiral 
speed induced by the temperature gradients: 



V = p'{~KdxT) 



13' sh _ h 



PsSdxTa{l - p' 



hM 
a 2e7 



(15) 

Continuing this analogy between the energy currents and 
charge currents, we can generalize the applicability of the 
result in Eq. (flSl) to transverse Neel DW— under the 
assumption of constant temperature gradients and van- 
ishing charge currents. Just like the (3 term is important 
for the current-driven DW dynamics,— the new viscous (3' 
term is important for the thermally induced DW motion 
below the Walker breakdown. 

To conclude, we phenomenologically introduced new 
/3— like viscous term into the LLG equation for the en- 
ergy currents. We speculate on a possibility of creating 
heat flows by microwave-induced periodic magnetization 
dynamics which should result in effective cooling of some 
regions, in analogy to the Peltier effect. To support it, we 
considered the DM spiral texture subject to charge cur- 
rent and found that the texture-dynamics induced heat 
flow is proportional to the Peltier coefficient and the new 
viscous coupling constant (3' . Thus, the materials with 
large Peltier coefficient and large (3' should be suitable 
for the realization of the microwave cooling by magneti- 
zation texture dynamics. The effective Peltier/Seebeck 
coefficient as well as the conductivity can be modified 
and tuned by the texture dynamics. These effects should 
be measurable in magnetic textures characterized by the 
length as small as ~ 10 nm and large viscous damping (3. 
Even in situations of pinned textures, the effects of EMF 
induced by viscous like term should be seen in mea- 
surements of the ac conductivity. We also conclude that 
the new viscous (3' term is important for the thermally 
induced DW motion. Bauer et alM- worked out very sim- 
ilar ideas for cooling by DW motion and thermoelectric 
excitation of magnetization dynamics. 
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